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We present results for the measurement of <j> meson production via its charged kaon decay channel 
<f) — > K + K~ in Au+Au collisions at ^/s NN = 62.4, 130, and 200 GeV, and in p + p and d+Au 
collisions at ^/s NN = 200 GeV from the STAR experiment at the BNL Relativistic Heavy Ion 
Collider (RHIC). The midrapidity (\y\ < 0.5) <j) meson transverse momentum (pr) spectra in central 
Au+Au collisions are found to be well described by a single exponential distribution. On the other 
hand, the pr spectra from p + p, d+Au and peripheral Au+Au collisions show power-law tails at 
intermediate and high pr and are described better by Levy distributions. The constant 4>/K~ yield 
ratio vs beam species, collision centrality and colliding energy is in contradiction with expectations 
from models having kaon coalescence as the dominant mechanism for <f> production at RHIC. The 
Q/<j> yield ratio as a function of pr is consistent with a model based on the recombination of thermal 
s quarks up to pr ~ 4 GeV/c, but disagrees at higher transverse momenta. The measured nuclear 
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modification factor, RdAu, for the <j> meson increases above unity at intermediate pr, similar to 
that for pions and protons, while Raa is suppressed due to the energy loss effect in central Au+Au 
collisions. Number of constituent quark scaling of both R cp and V2 for the <f> meson with respect 
to other hadrons in Au+Au collisions at % /s JVJV =200 GeV at intermediate pr is observed. These 
observations support quark coalescence as being the dominant mechanism of hadronization in the 
intermediate pr region at RHIC. 

PACS numbers: 25.75.Dw 



I. INTRODUCTION 

The (/>(1020) vector meson's properties and its trans- 
port in the nuclear medium have been of interest since 
its discovery The proper lifetime of the <ft meson is 
about 45 fm/c and it decays into charged kaons K + K~ 
with a branching ratio of 49.2%, and more rarely into the 
dilepton pairs e+e~ (B. R. of 2.97 x 1(T 4 ) and /Lt+ [T (B. 
R. of 2.86 x 1(T 4 ). 

The mechanism for <fi meson production in high en- 
ergy collisions has remained an open issue. As the 
lightest bound state of strange quarks (ss) with hidden 
strangeness, <j) meson production is suppressed in elemen- 
tary collisions because of the Okubo-Zweig-Iizuka (OZI) 
rule 0, 0, IH ■ The OZI rule states that processes with dis- 
connected quark lines in the initial and final state are sup- 
pressed. In an environment with many strange quarks, <j> 
mesons can be produced readily through coalescence, by- 
passing the OZI rule [f|. The <j> meson has been predicted 
to be a probe of the quark-gluon plasma (QGP) formed 
in ultrarelativistic heavy- ion collisions [1, 0, H, 0, US HH • 

On the other hand, a naive interpretation of 4> meson 
enhancement in heavy-ion collisions would be that the 4> 
meson is produced via KK —> in the hadronic rescat- 
tering stage. Models that include hadronic rescatterings 
such as RQMD and UrQMD [l]| have predicted an 
increase of the 4> to K~ production ratio at midrapid- 
ity as a function of the number of participant nucle- 
oids. This prediction was disproved in Au+Au collisions 

at V 1 ^ = 200 GcV ' from STAR y° ar 2001 data D3- 

With the higher statistics data newly recorded from the 
Solenoidal Tracker at RHIC (STAR) experiment [Hj], a 
precise measurement of the 4>/K~ ratio as a function of 
beam energy and collision centrality is presented to con- 
firm this finding in this paper. 

The in-medium properties of vector mesons in the hot 
and dense environment are also interesting [Ti|. The 
mass and width of the </> meson were predicted to change 
because of the partial restoration of chiral symmetry in 
the nuclear medium. Asakawa and Ko and Song [l~7| 
predicted that the <fi mass decreases as a result of many- 
body effects in a hadronic medium. A double <p peak 
structure in the dilepton invariant mass spectrum from 
relativistic heavy-ion collisions was proposed as a sig- 
nature of a phase transition from the QGP to hadronic 
matter [l8| . Other calculations have predicted that the <t> 
meson width can be widened significantly due to nuclear 
medium effects [ljl [2(| HJ- Recently, an interesting 4> 
mass modification at normal nuclear density in 12 GeV 
p + A interactions was observed in the dilepton channel 



(e+ + e") from the KEK experiment [HHH]. In STAR, 
we measure the <fi meson mass and width to compare with 
these predictions/measurements, using the decay channel 
4>^K + K-. 

From phcnomcnological analysis, it is suggested that 
the (j) meson mean free path in hadronic media is large 
because of its small cross section of scattering with 
hadrons [f|. Many other calculations also indicate that 
the <p meson has small rescattering cross sections with 
hadronic matter (l9l.[20|. However, after including three- 
and four-vector meson vertices into their hidden local 
symmetry model, Alvarez-Ruso and Koch [24|| found that 
the 4> meson mean free path in nuclear media is smaller 
than that usually estimated. Ishikawa et al. {25| pre- 
sented new data on near-threshold <p photoproduction on 
several nuclear targets. They found that the cross sec- 
tion between <fi meson and a nucleon, er^jv, is equal to 
35l^ mb, which appears to be much larger than previ- 
ous expectations although the experimental uncertainty 
is very large (26|. Meanwhile, Sibirtsev et al. (2?| pre- 
sented a new analysis of existing cf> photoproduction data 
and found g^n ~ 10 mb. Thus, u^n in heavy- ion colli- 
sions is still unclear. 

The measurement of collective radial flow (represented 
by {pt)) probes the equation of state of matter produced 
in nuclear collisions [H, [2j| . Strong radial flow has been 
observed for many particles such as 7r, K, and p(p) [30| . 
If the (f> meson has indeed small hadronic rescattering 
cross sections and decouples early from the collision sys- 
tem, contributions to the radial flow of the 4> meson will 
be mostly from the partonic stage instead of the hadronic 
stage. Thus the 4> meson may have a significantly smaller 
radial flow than other hadrons with similar mass, such 
as the proton, especially in central heavy- ion collisions. 
Therefore, (f> mesons may carry information about the 
conditions of nuclear collisions before chemical freeze-out. 
Thus it is important to experimentally measure and com- 
pare the freeze-out properties of the (j) to other hadrons 
as a function of centrality and collision species. A com- 
prehensive set of measurements will shed light on the 
characteristics of 4> meson production and the evolution 
of the collision system. 

The elliptic flow parameter v 2 is a good tool for study- 
ing the system formed in the early stages of high energy 
collisions at RHIC. It has been found that at low px 
(0 < pr < 2 GeV/c), the dependence of «2 on particle 
mass |3j HU, HH is consistent with hydrodynamic cal- 
culations in which local thermal equilibrium of partons 
is assumed [34|, HH, H(| H3, Hl|. This observation indi- 
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cates that thermally equilibrated partonic matter may 
have been created at RHIC. However, at intermediate 
Pt (2 < pr < 5 GeV/c), the measured v 2 for various 
hadrons seems to depend on the number of constituent 
quarks in the hadron rather than its mass, consistent 
with the results from coalescence/recombination mod- 
els [H El EH, E2, S3, El- Since 4> is a meson but has a 
mass close to that of the proton and A, the measurement 
of the 4> meson elliptic flow will provide a unique tool for 
testing the above statement. 

Current measurements of various hadrons by STAR 
(A, p, A'g, #(892)*, h ± , etc.) show that the nuclear 
modification factor R cp for baryons differs from that of 
mesons [5^, EH , consistent with the prediction of quark 
coalescence/recombination models [H, EE E3, EH- Be- 
cause of the value of the (f> meson mass, a comparison of 
R cp for the <fi with these previous measurements will con- 
clusively determine if the observed difference is driven by 
particle mass or particle type. 

In previous studies, the production of the meson 
has been measured in Au+Au collisions at ^/s Niv = 
130 GeV [Ii| and in Au+Au and p + p collisions at 
v /s^7 = 200 GeV [3 at RHIC. The Au+Au 200 GeV 
data presented in this paper were taken from the year 
2004 run, where the number of events is approximately 
ten times larger than the previously reported number for 
the Au+Au run in year 2001 [l4j|. The data from the 
year 2004 run contain the data reported in ref . [H(| [Hl| ■ 
It has been found that the results from the two runs are 
consistent with each other. In this paper, we present 
systematic measurements of 4> meson production over a 
broad range of collision energies and system sizes, in- 
cluding Au+Au collisions at t/s~^ = 62.4, 130 [49l |. 
and 200 GeV, and p + p [3] and d+Au collisions at 
Y / s^" = 200 GeV from the STAR experiment. In Sec. 
II, we briefly introduce the STAR detector and discuss 
our analysis method (event-mixing technique) in detail. 
In Sec. Ill, we present the measurement of </> meson in- 
variant mass distributions (III A), transverse mass m<r) 
spectra (III B, particle ratios (III C), nuclear modifica- 
tion factors (III D) and the elliptic flow parameter v 2 
(III E); we also discuss the physics implication of each of 
these results. A summary and conclusions are presented 
in Sec. IV. 

II. DATA ANALYSIS 

A. Experimental Setup 

The STAR detector [IH[ consists of several subsystems 
in a large solcnoidal analyzing magnet. We discuss here 
the subdetectors used in the analyses relevant to this pa- 
per. With its axis aligned along the beam direction, the 
time projection chamber (TPC) [52| is the main tracking 
device for charged particles, covering a pseudorapidity 
range |?7| < 1.8 and providing complete azimuthal cover- 
age. The entire TPC is located inside a solcnoidal mag- 
net, and data are taken at the maximum magnetic field 
\B Z \ = 0.5 Tcsla, where the z axis is parallel to the beam 



direction. Radial-drift TPCs (FTPCs) [H are also in- 
stalled to extend particle tracking into the forward and 
backward regions (2.5 < < 4.0). Surrounding the 
TPC is the central trigger barrel (CTB) [HI], which is a 
scintillator counter array whose analog signal is sensitive 
to the total charged particle multiplicity with coverage 
|?7| < 1.0. A pair of beam-beam counters (BBCs) at 
3.3 < n < 5.0 and a pair of zero degree calorimeters 
(ZDCs) (HHJ at 9 < 2 mrad are located on either side 
of the collision region along the beam line, and are used 
to provide event triggers for data taking. A more de- 
tailed description of the STAR detector can be found in 
Ref. [TH| and references therein. 

B. Event selection 

1 . Trigger selection 

The Au+Au data used in this analysis were taken with 
two different trigger conditions. One was a minimum-bias 
(MB) trigger requiring only a coincidence between both 
ZDCs. The other was a central trigger additionally re- 
quiring both a large analog signal in the CTB indicating 
a high charged particle multiplicity at midrapidity and 
a small ZDC signal. The ZDCs measure beam-velocity 
neutrons from the fragmentation of colliding nuclei and 
were used as the experimental level-0 trigger for select- 
ing d+Au and Au+Au collisions for their respective runs. 
For p + p data taking, the BBCs were used as trigger 
detectors. The central trigger corresponds to approxi- 
mately the top 15% and 12% of the measured cross sec- 
tion for Au+Au collisions at 130 GeV and 200 GeV, re- 
spectively. Data from both the MB and central triggers 
were used for this analysis. For the Au+Au 62.4 GeV 
data set, only MB triggered events were used. For p + p 
collisions at 200 GeV, the MB trigger was used in the 
analysis. It was based on a coincidence between the two 
BBCs. The BBCs are sensitive only to the non-single 
diffractive (NSD) part (30 mb) of the p+p total inelastic 
cross section (42 mb) [Ffj . 

1 . Vertex cuts and centrality selection 

The longitudinal z position of the interaction point is 
determined on-line by the measured time difference of the 
two ZDCs's signals. A cut on the z position of the in- 
teraction point is applied on-line for all data sets (except 
p + p) in order to maximize the amount of useful data for 
physics analysis, since events with primary vertices far 
away from the center of the TPC have a significant non- 
uniform acceptance. In the off-line data analysis further 
cuts are applied on the z position of the reconstructed 
primary vertex (Vz), to ensure nearly uniform detector 
acceptance. These cuts are listed in Table |U 

To define the collision centrality for the Au+Au data, 
the raw charged hadron multiplicity distribution in the 
TPC within a pseudo-rapidity window < 0.5 (\rj\ < 
0.75 was used for the Au+Au 130 GeV data) was di- 
vided into several bins. Each bin corresponds to a cer- 
tain fraction of the total inelastic cross section [56[ . For 
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System 


\/ S NN 

(GeV) 


Trigger 


\Vz\ 
(cm) 


Centrality 


Events 


Year 


Au+Au 


62.4 


MB 


< 30 


0-80% 


6.2 x 10 6 


2004 


Au+Au 


130 


MB 


< 80 


0-85% 


7.6 x 10 5 


2000 


Au+Au 


130 


Central 


< 80 


0-11% 


8.8 x 10 5 


2000 


Au+Au 


200 


MB 


< 30 


0-80% 


1.4 x 10 7 


2004 


Au+Au 


200 


Central 


< 30 


0-5% 


4.8 x 10 6 


2004 


P + P 


200 


MB 


< 50 


MB (NSD) 


6.5 x 10 6 


2002 


d+Au 


200 


MB 


< 50 


MB 


1.4 x 10 7 


2003 



TABLE I: Data sets used in the analysis. Cuts on Vz, the 
selected centrality ranges and the final number of events in- 
cluded in the analysis after all cuts/selections are also shown. 



the d+Au data, the raw charged hadron multiplicity in 
the east (Au-direction) FTPC (-3.8 < ?y < -2.8) was used 
for the centrality definition to avoid auto-correlation be- 
tween centrality and the measurements of charged par- 
ticles at midrapidity in the TPC [5(|. We defined four 
centrality bins for the Au+Au 62.4 GeV data (0-20%, 
20-40%, 40-60%, 60-80%), three centrality bins for the 
Au+Au 130 GeV data (0-11%, 11-26%, 26-85%), nine 
centrality bins for the Au+Au 200 GeV data (0-5%, 0- 
10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, 
70-80%) and three centrality bins for the d+Au 200 GeV 
data (0-20%, 20-40%, 40-100%). The 80-100% most pe- 
ripheral Au+Au collision data were not used because of 
the rapidly decreasing trigger and vertex finding efficien- 
cies for low multiplicity events. Table U lists the data sets 
used along with centrality selections and final numbers 
of events after these cuts. Note that the elliptic flow fa) 
measurement only from 200 GeV MB Au+Au data is pre- 
sented in this paper, as the statistics for the i>2 analysis 
is not sufficient for the 62.4 and 130 GeV Au+Au data. 

C. Track selection and particle identification 

1 . Track selection 

Several quality cuts were applied to ensure selection 
of good tracks. During the TPC track reconstruction, a 
charged track was extrapolated back to the beam line by 
using the reconstructed helix parameters. If the distance 
of closest approach (DCA) of the track to the event ver- 
tex was less than 3 cm and the track had at least ten hit 
points in the TPC, the reconstructed track was labeled 
as a primary track. The helix parameters for primary 
tracks were then refitted by requiring that the helix pass 
through the primary vertex location. This procedure im- 
proved the momentum resolution of tracks. Since the 
<p meson has a very short lifetime, it decays at the pri- 
mary vertex position. Thus only primary tracks were 
used for the meson analysis. As a systematic check, 
the DCA selection for primary tracks was changed from 
3 cm to 1 cm. The differences in the results were small 
and were included in the estimate of systematic uncer- 
tainties. Tracks with transverse momentum less than 0.1 




momentum x charge (Gev/c) 



FIG. 1: (Color online) Measured (dE/dx) vs momentum X charge 
of reconstructed tracks in the TPC. The figure is generated from 
Au+Au 62.4 GeV data. 



GeV/c were not used, as their combined acceptance and 
efficiency becomes very small. Each track included in the 
<j> analysis was required to have at least 15 hit points out 
of 45 used in the fitting of the tracks helix parameters. 
The ratio of the number of space points used in the track 
reconstruction to the maximum possible number of hit 
points was required to be greater than 55% to avoid split 
tracks where a real track is reconstructed in two or more 
segments. A pseudo-rapidity cut |^| < 1.0 (\n\ < 1.1 for 
Au+Au 130 GeV data) was applied to select tracks that 
are well within the TPC acceptance. 

2 . Kaon selection 

Particle identification (PID) was achieved by correlat- 
ing the ionization energy loss {dE/dx) of charged par- 
ticles in the TPC gas with their measured momentum. 
The measurement of mean dE/dx was achieved by aver- 
aging the measured dE/dx samples along the track after 
truncating the top 30%. The measured (dE/dx) ver- 
sus momentum curve is reasonably well described by the 
Bethe-Bloch function [57} smeared with the detector's 
resolution (note that the Bichsel function was used to fit 
the {dE/dx) plot in Au+Au 200 GeV from the year 2004 
run 58]). The n„ values for kaons are calculated via 



I measured 



1 (dE/dx) 

R dE I dx eX pectea 



(i) 



where (dE / dx) measured and dE / dx expect ed are (dE/dx) 
measured by TPC and calculated analytically, respec- 
tively, and R denotes the dE/dx resolution of the track 
which is found to range between 6% and 10%. R is de- 
termined experimentally and depends on the event multi- 
plicity and the number of dE/dx samples from the track 
used to calculate the mean value. Tracks within 2a of 
the kaon Bcthc-Bloch curve were selected as kaon can- 
didates. Figure [1] presents the measured (dE/dx) versus 
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momentum x charge in Au+Au collisions at 62.4 GeV. 
Table HT1 lists all the track cuts applied in the analysis. 



Cut parameter 



Value 



Track DCA (cm) 

Track N Flt 

Track N F it/N Max 

Track momentum (GeV/c) 

Track transverse momentum (GeV/c) 

Kaon dE/dx 

4> candidate's S— dip— angle (radians) 
(j> candidate's rapidity 



< 3 
> 15 

> 0.55 
0.1 < p < 10 
0.1 < p T < 10 

|ricr| < 2.0 (for kaon) 

> 0.04 

\y\ < 0.5 (for spectra) 
\y\ < 1.0 (for v 2 ) 



TABLE II: Track cuts used in the analysis, where Nfu and 
Nmux represent the number of fitted hits and the maximum 
number of hits for TPC tracks, respectively. 

Note that from the (dE/dx) measurement, kaons can- 
not be clearly separated from pions above p ^0.6 GeV/c 
and from protons/antiprotons above p ~ 1.1 GeV/c. Also 
note that the electron and positron dE/dx bands cross 
the bands for pions, kaons, and protons/antiprotons. 
Therefore selected kaon candidates are contaminated by 
electrons/positrons, pions and protons/antiprotons vary- 
ing with p. Contamination by these charged particles 
in the kaon sample brings in additional real correlations 
(such as particle decays) which cannot be subtracted by 
the event-mixing method. We have varied the n a cut 
for kaon to investigate the efficiency and combinatorial 
background dependence of the <p signal extraction on this 
cut. The resulting systematic uncertainties have been in- 
cluded in the estimate of the total systematic errors. 

D. Event mixing and raw yield extraction 

The meson signal was generated by pairing all 
K + K~ tracks from the same event that passed the selec- 
tion criteria and by then calculating the invariant mass 
rriiny for all possible K + K~ pairs. As random combina- 
tions of K + K~ pairs are dominant in this process, the 
resulting same-event invariant mass distribution contains 
the 4> meson signal on top of a large combinatorial back- 
ground. An event-mixing technique [H, [6(| was applied 
to calculate the shape of the combinatorial background, 
where the invariant mass was calculated by pairing two 
kaons from two different events with same primary vertex 
and multiplicity bins (mixed event). Ideally, since it com- 
bines two different events, the mixed-event distribution 
contains everything except the real same-event correla- 
tions. 

The STAR TPC has symmetric coverage about the 
center of the collision region. However, variations in the 
acceptance occur, since the collision vertex position may 
change considerably event-by-event. This variation in the 
collision vertex position gives rise to a nonstatistical vari- 
ation in the single-particle phase-space acceptance, which 



would lead to a mismatch between the mixed-event and 
same-event invariant mass distributions. This mismatch 
would prevent the proper extraction of the <f> meson sig- 
nal. By sorting events according to their primary ver- 
tex Vz position and performing event-mixing only among 
events in the same vertex bin, the mismatching effect is 
minimized. In this analysis, events were divided into Vz 
bins that were 6 cm wide (Vz resolution is ~ 0.3 mm) 
for event-mixing. To further improve the description of 
the background, two events were only mixed if they had 
similar event multiplicities. These requirements ensure 
that the two events used in mixing have similar event 
structures, so the mixed-event invariant mass distribu- 
tion can better represent the combinatorial background 
in the same-event invariant mass distribution. Consis- 
tent results were obtained when we constructed the back- 
ground distribution using like-sign pairs from the same 
event. 

To reduce statistical uncertainty in the mixed event, 
each event was mixed with five to ten other events (de- 
pending on the collision system) . To extract the (j> meson 
signal, the mixed-event and same-event K + K~ invari- 
ant mass distributions were accumulated and the mixed- 
event distribution normalized to the same-event distribu- 
tion in the region above the 4> mass, 1.04 < rrii nv < 1.06 
GeV/c 2 , and subtracted in each px and y (rapidity) bin 
for every collision centrality. We varied the normalization 
region and normalization factor to estimate the system- 
atic uncertainty on the normalization, and the estimated 
uncertainty was included in the quoted total systematic 
uncertainty. 

Despite the requirements for mixing events described 
above, a residual background remains over a broad mass 
region in the subtracted invariant mass distribution. This 
is due to an imperfect description of the combinatorial 
background and the fact that the mixed event cannot 
account for the real correlated background from decay 
pairs due to Coulomb interactions, photon conversions 
(7- 



~e ), and particle decays such as K a * — > K + tt 



p u — > 7r + 7r~, Kg — > 7r + 7r~, and A — ► pw~ [61| . For exam- 
ple, when both pions from a K s decay are misidentified 
as kaons, the real correlation from decay will remain in 
the same-event as a broad distribution, but will not be 
reproduced by the event-mixing method. 

Due to contamination of electrons / positrons in the se- 
lected kaon sample, the K + K~ invariant mass distri- 
bution contains residual background near the threshold 
from correlated e + e _ pairs, mainly from photon conver- 
sions (7 — > e + + e~). The 5-dip-anglc between the photon 
converted electron and positron is usually quite small. 
The (5-dip-angle is calculated from 



8— dip— angle 



-Ir gTlPTg + PzlPz2 , 

P1P2 



(2) 



where pi,P2, Pti, PT2, Pzi> Pz2 are momentum and trans- 
verse and longitudinal momentum components of the two 
tracks; this represents the opening angle of a pair in the 
Pz~Pt plane. We required the 5-dip-angle to be greater 
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FIG. 2: (Color online) Background-subtracted invariant mass dis- 
tributions at 0.4 < p T < 0.8 GeV/c in d+Au 200 GeV collisions 
(0-100%) with (solid points) and without (open points) the (5-dip- 
angle cut. The dashed curves show a Breit-Wigner (see the text 
for details) + linear background function fit to the case with the 
<5-dip-angle cut. 



than 0.04 radians. This cut was found to be very effective 
in removing the photon conversion background while only 
reducing the cp reconstruction efficiency by ~ 12%. Fig- 
ure [2] shows two background-subtracted invariant mass 
distributions with and without the (5-dip-angle cut. One 
can see that the peak from photon conversion(mi n/u < 1.0 
GeV/c 2 ) is very effectively removed by this cut. 

Figure [3] shows the K + K~ invariant mass distribu- 
tions for p+p collisions at 200 GeV [(a) and (b)], 60-80% 
Au+Au collisions at 62.4 GeV [(c) and (d)] and 0-10% 
Au+Au collisions at 200 GeV [(e) and (f)]. Solid circles 
in the upper panels are same-event pairs, whereas the his- 
tograms are from mixed-event pairs. The cp meson peak 
is clearly visible for p+p 200 GeV and Au+Au 62.4 GeV 
(60-80%) in Figs.[3ja) and (c) before background subtrac- 
tion, but not for Au+Au 200 GeV (0-10%) [Fig. He)] due 
to its smaller signal significance. However, after back- 
ground subtraction, the cp mass peak can be seen clearly 
for all data sets. The lower panels in Fig. [3] show the 
mixed-event background-subtracted cp invariant mass dis- 
tributions. Raw yields for the cp meson were determined 
by fitting the background-subtracted n%i nv distribution 
with a Breit-Wigner function superimposed on a linear 
(or polynomial) background function 



dN 



AT 

—^- f ^+B(rn inv ), (3) 



where A is the area under the peak corresponding to the 
number of cp mesons, T is the full width at half maximum 
(FWHM) of the peak, and mo is the resonance mass posi- 
tion. B(rrii nv ) denotes a linear [B(rrii nv ) = po + p\mi nv , 
shown by a dashed line in Figs. GUb), (d) and (f)] or 
polynomial [B(m inv ) = p Q + pim lnv + pim 2 inv , shown by 



a dot-dashed line in Fig. [3jf )] residual background func- 
tion. The parameters po, p\ and pi of B{m inv ) and A, 
mo and T are free parameters. 

E. Efficiency correction 

The cp acceptance and reconstruction efficiency were 
calculated using an embedding technique, in which sim- 
ulated tracks were embedded into real events. The num- 
ber of embedded simulated tracks is approximately 5% 
of the measured multiplicity of the real event. The cp me- 
son decay (cp — ► K + K~) and the detector responses were 
simulated by the GEANT program package [62j and the 
simulated output signals were embedded into real events 
before being processed by the standard STAR event re- 
construction code. Embedded data were then analyzed 
to calculate tracking efficiencies and detector acceptance 
by dividing the number of reconstructed cp by the num- 
ber of input cp in the desired kinematic regions. Figure 0] 
shows examples of correction factors (tracking efficiency 
x acceptance) for our analysis as a function of cp me- 
son pt for selected centrality bins for Au+Au 200 GeV, 
d+Au 200 GeV and p + p 200 GeV collisions. It can be 
seen that the overall correction factors increase from a 
few percent at low pr to over 30% at high px- haw ef- 
ficiency at low pt is mainly due to poor acceptance of 
the daughter tracks. The efficiency is lower in more cen- 
tral collisions because of the increasing occupancy in the 
TPC [H. 

F. Vertex finding and trigger efficiency 
correction 

For p + p and d+Au data, the trigger efficiency is less 
than 100% [63j. The MB trigger for p + p data was found 
to trigger ~ 87% of p + p NSD events. For d+Au data, 
the trigger efficiency was found to be ~ 95%. These 
trigger efficiencies were used to normalize the measured 
yield and the corresponding uncertainties are added to 
the total systematic errors for p+p and d+Au data. The 
MB trigger efficiency for Au+Au data is essentially 100% 

It was found that the event vertex finding efficiency, 
which is the fraction of events having reconstructed ver- 
tices, drops rapidly for low multiplicity events [5(J. For 
d+Au collisions, the vertex efficiency was 88% for the 
most peripheral bin (40-100%), 93% for the MB events 
and 100% for the central and middle central bins (0-20% 
and 20-40%). For Au+Au collisions, the vertex finding 
efficiency was 99.9% and pr independent [64| due to the 
increased track multiplicity in those collisions. However 
the overall vertex finding efficiency was found to be 98.8% 
for the MB p + p data by applying an additional BBC 
time difference selection (i.e., the NSD requirement), and 
the effect of correction was negligible in this analysis. 

G. f2 measurement 

1 . Reaction plane method 
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FIG. 3: (Color online) Upper panels: same-event (full points) and mixed-event (solid line) K + K invariant mass distributions at 
0.6 < p T < 1.4 GeV/c in p + p 200 GeV collisions (a), 0.8 < p T < 1.2 GeV/c in Au+Au 62.4 GeV collisions (60-80%) (c) and 
0.8 < pt < 1.2 GeV/c in Au+Au 200 GeV collisions (0-10%) (e). Lower panels: the corresponding <j> meson mass peaks after subtracting 
the background. Dashed curves show a Breit-Wigner + linear background function fit in (b), (d). In (f), both linear and quadratic 
backgrounds arc shown as dashed and dot-dashed lines, respectively. 
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FIG. 4: (Color online) Reconstruction efficiency including accep- 
tance of eft meson as a function of Pt in several centrality bins of 
Au+Au, d+Au and p + p 200 GeV collisions. 



E 



d 3 N 



1 d 2 N 



d 3 p 2ir ptdptdy 



2v n cos[?i(<^ — \& r )] 



(4) 

where ^ r is the real reaction plane angle and tp is the par- 
ticle's azimuthal angle. The coefficient V2 in the second- 
order term of the expansion is the dominant part and is 
called the second harmonic anisotropic flow parameter, 
or elliptic flow. 

The real reaction plane angle is not known, but 
can be estimated experimentally (67j . In our analysis, 
the estimated reaction plane angle from the second-order 
harmonic (^2) was used. This has a finite resolution 
due to a limited number of particles available in each 
event and a different event-by-event V2 , which is used for 
the estimation. The estimated reaction plane resolution 
was used to correct the observed to obtain the final 
estimation of i>2. 

The event plane angle "J" 2 was calculated by the equa- 
tion 



^2 



■ tan 



J2j WjSm.(2(pi) 
J2i Wi cos(2c/Pi) 



(5) 



We employed the STAR standard reaction plane 
method as described in Refs. [g^ |66| . which uses a 
Fourier expansion to describe particle emission with re- 



where the sums are over all charged particles used for 
reaction plane determination, Wi and (fi are the weight 
and azimuthal angle for the ith particle in a given event, 
respectively. The weights include both a pt weight and 
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if weight. The px weight was taken to be the particle 
Pt up to 2.0 GeV/c and constant (2.0) above that [6(| . 
The f weight was taken to be the reciprocal of the ip 
distribution (normalized by the average entries) for all 
selected tracks. The autocorrelations were eliminated by 
excluding all kaon tracks used in the </> invariant mass 
calculation from the reaction plane angle estimation . 
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FIG. 5: Event plane ^2 resolution as a function of centrality in 
Au+Au 200 GeV collisions, where the vertical axis starts from 0.3 
for clarity. 

The reaction plane resolution was then calculated by 

(cos[2(<F 2 - 9 r )]) = C(cos[2(^ - * r )]), (6) 

where "F^ is the calculated reaction plane angle of the 
subcvcnt, and C is a constant calculated from the known 
multiplicity dependence of the resolution [66| . This reso- 
lution was determined by dividing each event into two 
random subevents, a and b, with equal multiplicities. 
The reaction plane resolution therefore corresponds to 
how accurately the event plane angle represents the real 
reaction plane; due to its the definition in Eq. ([6]), a value 
of unity indicates ideal resolution. Figure [5] presents the 
reaction plane resolutions in different centrality bins for 
Au+Au 200 GeV collisions. 

The combinatorial background in the <f> invariant mass 
distribution was also calculated by an event-mixing tech- 
nique as described above. To guarantee the mixed-event 
sample would represent the combinatorial background, 
an additional cut, the reaction plane angle difference be- 
tween two mixed events, was required to be less than 
0.1-7T rad in the event-mixing procedure. After back- 
ground subtraction, <f> meson yield was extracted in each 
{pT, f ~ ^2) bin. The yield distribution as a function of 
ip — ^2 was fitted by the function 



(7) 



A[ 1 + 2vf s cos[2(v? - # 2 )] 



result for the Au- 
Fig.l] 



Au data at 200 GeV is presented in 
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FIG. 6: ip — 3?2 distribution for cf> meson at 1.5 < pt < 2.0 GeV/c 
in Au+Au collisions (0-80%) at 200 GeV. The line is the fitting 
result. Error bars are statistical only. 

The measured v^" 3 was then divided by the reaction 
plane resolution to obtain the final V2, i.e., 



,obs 



l'2 



(coB[2(tt a - *r)]> 



(8) 



to extract the v^ s value, where A is a constant. A typical 



Simulation studies have found that the measured vi is 
about 7% (relative to the real f 2 ) lower than the real V2 
due to binning effects (five bins in if — VE^); a correction 
has been made to the measured V2 to account for this 
effect. 

2 . Invariant mass method 

A new method, namely, the invariant mass method. 
was also used to extract the elliptic flow v-i of the <p me- 
son. The method was proposed in Ref. (68|, which de- 
composes the anisotropic flow v n of a short-lived particle 
from that of all possible daughter pairs as a function of 
invariant mass. For extracting the vi of the </> meson, it 
utilizes the fact that the t>2 of K + K~ pairs is composed 
of the V2 of the combinatorial background and the vi 
of the <f> meson. Following the mixed-event background- 
subtraction procedure described in Sec. II D, the number 
of K + K~ pairs in each invariant-mass bin were counted, 
irrespective of the pair azimuth. Then 

in v 

) + N B (m inv ), (9) 

where and Nb are from the <j> signal and the back- 
ground, respectively. Once has been extracted via 
event-mixing and fitting the <f> mass peak with Eq. (fS]) 
for each pt bin as discussed in Section II D, Nb can be 
obtained from Eq. ©. 
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The same-event v 2 for K + K pairs vs. 
can be described by the function 



invariant mass 



) = a(m inv )v2S + [1 - a(rninv)]v2B(m inv ), (10) 

where ^(^nn) is the v 2 of same-event K + K~ pairs, 
V 2S = v 2<t> is the v-x of the meson, v 2 b is the effec- 
tive V2 of the combinatorial background and a(rrii nv ) = 
N^{mi nv )/N K + K -(mi nv ) is the ratio of the signal to 
the sum of the background and (f> signal. The reaction 
plane angle \E , 2 was estimated in the same way as for the 
reaction plane method described in the previous section. 
Therefore the two methods arc not completely indepen- 
dent. V2(rrii nv ) can then be calculated from the following 
equation (66| for each m inv bin 



v 2 (m inv ) = (cos[2(<p 



KK 



*a)]>, 



(11) 



where (Pkk is the azimuthal angle of the K + K~ pair. 

Under the assumption that the background contribu- 
tion to V2(jrii nv ) [the second part on right side of equa- 
tion pop ] is smooth as a function of rrii nv [68j . a poly- 
nomial function, po+pimi nv +p2irq nv , can be used to 
parametrize the background V2B vs rrii nv . V2S is then ob- 
tained by fitting v 2 by Eq. (|10[) in each px bin, with V2S as 
a free parameter. Figure [7] shows (cos[2(ip K + K - — ^2)]) 
[i.e. V2 in Equation[TT] vs m inv for 0.5 < pr < 1.0 GeV/c 
in Au+Au 200 GeV collisions (0-80%), where the solid 
curve is the result of fitting Eq. (flU)) . At the same time, 
(sm[2(tpfc+K- — ^2)]) vs rrii nv (open points) is found, as 
expected, to be consistent with zero due to collisional ge- 
ometry symmetry [§6|. The V2S value [i.e., v% bs in Eq. J7J] 
determined by the fit was corrected for the reaction plane 
resolution to obtain the final V2 for the <p meson. The fi- 
nal V2 results and related discussions will be presented in 
Sec. Ill E. 

H. Systematic uncertainties 

Major contributions to the systematic uncertainties 
come from variations in the procedure used for extract- 
ing the yields from the K + K~ invariant mass distribu- 
tions and from variations in the determination of tracking 
and particle identification efficiencies. Different residual 
background functions (first-order vs. third-order poly- 
nomial curves) were used to estimate the uncertainty of 
the raw yield extraction in each bin, and it was found to 
be of the order of ~ 4.5%. The uncertainty due to dif- 
ferent mixed-event normalization factors was estimated 
to be ~ 2.1% by varying the normalization region in the 
mixed-event background distribution. The uncertainty 
from tracking and PID efficiencies was estimated to be 
~ 8%, by varying the kinematic and PID cuts on the 
daughter tracks. 

The overall systematic uncertainty was estimated to 
be approximately 10% for the yield (dN/dy), and 10% 
for (px) for the Au+Au and the d+Au data. It includes 
an additional contribution from the difference between 
exponential and Levy function fittings of the transverse 
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FIG. 7: (Color online) (cos[2(ip K+ K - — ^2)]) (full rc d points) and 
(sin[2(ipj r + K - — ^2)]) (open blue points) as a function of mi nv 
of K+K~ pairs at 0.5 < p T < 1.0 GeV/c in Au+Au 200 GeV 
collisions (0-80%), where the solid curve is the result of fitting by 
Equation (1 1 1 . The arrow shows the position of <f> invariant mass 
peak. The dash line shows zero horizontal line. 



mass or transverse momentum distributions. The sys- 
tematic uncertainty in the overall normalization for the 
p + p 200 GeV data was found to be 15% for dN/dy and 
5% for (pt)j including uncertainties due to vertex finding 
and trigger inefficiency for low multiplicity events. 

Systematic uncertainties for the 1*2 measurement from 
the two different v 2 extraction methods show px and cen- 
trality dependences, which mainly result from the de- 
termination of S/(S+B) ratios for the invariant mass 
method and from the removal of residual background in 
the reaction plane method, respectively. In our analysis, 
the point-to-point systematic errors included contribu- 
tions from the following: 

(i) Difference in finding the 0-meson signal via bin-by- 
bin counting or Breit-Wigncr function fitting meth- 
ods; 

(ii) Difference due to the residual background fitting 
function: first- or third-order polynomial functions; 

(iii) Difference in combinatorial background determina- 
tion: rotation of the background (the mixed event 
is from the azimuthal angle rotation of all tracks 
from the same event) or event-mixing (the current 
method); 

(iv) Difference in v 2 calculation: centrality-by- 
centrality V2 calculation and then weighting to 
get the final MB v 2 or direct calculation of the V2 
through MB raw yield fitting. 

III. RESULTS 

A. Mass and width 
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FIG. 8: (Color online) Masses and widths (FWHMs) of <j> as a function of p T in p + p 200 GeV (NSD), ci+Au 200 GeV (0-20%), Au+Au 
62.4 GeV (0-20%) and Au+Au 200 GeV (0-5%) collisions, with the corresponding PDG values. 



Figure [8] shows the </> invariant mass peak position and 
width (FWHM) as a function of p T for Au+Au 200 GeV 
(0-5%), Au+Au 62.4 GeV (0-20%), d+Au 200 GeV (0- 
20%) and p + p 200 GeV (NSD) collisions. In the larger 
Pt region (> 1 GeV/c), the measured mass and width for 
the 4> meson are consistent with those from Monte Carlo 
(MC) embedding simulations in various collision systems 
and at different energies. At low pr (< 1 GeV/c), the 
measured </> meson mass is lower and the width is larger 
than from simulation. The drop of the (j) meson mass 
in both real data and simulation at low pt is due to 
the multiple scattering energy loss of low pt tracks in 
the detector, which is not fully corrected during track 
reconstruction. 

Figure [9] shows shape comparisons between experimen- 
tal and MC invariant mass distributions for the 4> meson 
at 0.6 < p T < 1.0 GeV/c in p+p 200 GeV (NSD) and 
Au+Au 200 GeV (0-5%) collisions. The real data in- 
variant mass peak (solid circles) is wider than that from 
standard MC data set (1) (open circles). If the momen- 
tum resolution for low pt kaons used in the simulations 
is increased by 50%, e.g., kaon momentum resolution at 
350 MeV/c increases from -2% [MC data set (1)] to -3% 
[MC data set (2)], the (f> meson width from simulation re- 
produces the measured width from real data as shown by 
open diamonds. This decreased momentum resolution in 
MC could be possible considering uncertainties in simula- 
tions for the amount of material between the TPC active 
volume and the primary collision vertex and residual ge- 
ometry alignment issues. These remaining issues for the 



differences in mass and width of </> mesons between real 
data and simulations have limited our sensitivity to pos- 
sible small modifications of </> meson properties in the 
medium produced at RHIC collisions. It should also be 
noted that to really trace down the possible modifica- 
tion of the 4> meson mass and width in heavy-ion colli- 
sions, measurements through the dilepton decay channel 
are needed. An interesting excess on the low-mass side 
of the <f> meson invariant mass peak was observed by an 
e + + e~ channel in the low f3j region ((3j < 1.25) for 
12 GeV jP+Cu interactions from the recent KEK exper- 
iment [22l [23[ . This may indicate a vector meson mass 
modification at normal nuclear density. <fi measurements 
using the dilepton channel will hopefully be addressed in 
STAR in year 2010 with the time-of-flight detector up- 
grade under construction. 

B. Spectra 

4> meson differential invariant yields were calculated 
by correcting the extracted raw yield by tracking ef- 
ficiency, detector acceptance and the decay branching 
ratio. Figure [TU] shows the <p meson transverse mass 
(jtit = y/ p T + tuq, where mo is the mass of <f> meson) dis- 
tributions from Au+Au collisions at y/s NN = 62.4, 130, 
and 200 GeV and from p + p (NSD) and ti+Au collisions 
at y / s^ N " = 200 GeV. All spectra are from midrapidity, 
\y\ < 0.5, with pt coverage above 0.4 GeV/c. For clarity, 
distributions for different ccntralitics are scaled by fac- 
tors indicated in the figure. Lines in the figure represent 
fits to the transverse mass distributions for different ccn- 
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FIG. 9: (Color online) Invariant mass distributions of <j> meson at 
0.6 < p T < 1-0 GeV/c in p + p 200 GeV (NSD) and Au+Au 
200 GeV (0-5%) collisions. Solid symbols: experimental data. 
Open symbols: MC simulation. Curves are the results of a Breit- 
Wigner function fit. Note: Two sets of MC data are shown for 
Au+Au 200 GeV, and see text for details. 



tralities. The 62.4, 130 and 200 GeV Au+Au data were 
fitted by the exponential function 



1 



d 2 N 



dN/dy 



-(mT—mo)/Te: 



2irni T dm T dy 2irT exp (m + T exp ) 

(12) 

where the slope parameter T exp and yield dN / dy are free 
parameters. For p+p 200 GeV and <i+Au 200 GeV data, 
the distributions were fitted by a Levy function [6^ . F70[ | 
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where n, the slope parameter TL evy , and yield dN/dy, are 
free parameters. For the four most peripheral ccntrality 
bins (40-50%, 50-60%, 60-70% and 70-80%) in Au+Au 
collisions at 200 GeV, the distributions were better fit by 
a Levy rather than by an exponential function. In fact, 
the exponential function [Eq. (fT2|) ] is the limit of the Levy 



function [Eq. (|13p ] as n approaches infinity; i.e., T ( 



^Levy {n—>oo). Table ITTTl lists the extracted slope parame- 
ter T, mean transverse momentum (pr), and yield dN/dy 



from the best fits to the spectra. Overall estimated sys- 
tematic uncertainties on these quantities are also listed. 

The <j) meson transverse mass spectra in central 
Au+Au collisions can be well described by a single 
m^-exponential function, while the spectra in ci+Au, 
p + p and peripheral Au+Au are better described by 
a Levy function, due to the power-law tail at interme- 
diate and high pt- Figure [IT] compares the transverse 
momentum spectra shapes in different 200 GeV colli- 
sion systems (0-5% Au+Au, 0-20% d+Au, and inelas- 
tic p + p). The spectra arc normalized by the num- 
ber of binary collisions (N i n ) and number of partici- 
pant pairs (N part /2). N oin and N part were determined 
by Glauber model calculations [5(|. We again point out 
that STAR only triggered on NSD p + p events (mea- 
sured onsd = 30.0 ± 3.5 mb) [7l|, while the Glauber 
model calculations use the p + p inelastic cross section 
{pinei —42+ 1.0 mb). Thus the NSD p + p spectrum was 
normalized to the inelastic yield by a correction factor of 
30/42. 

A change in the shape of spectra from p + p, d+Au 
and peripheral Au+Au collisions to central Au+Au col- 
lisions is clearly visible. In comparison with the fit- 
ting result for 200 GeV p + p collisions in the high pr 
(4.0 < pr < 6.0 GeV/c) region, the Nb in normalized yield 
is suppressed in central Au+Au collisions at 200 GeV, 
while no suppression is observed for d+Au collisions at 
200 GeV. Since particles with high transverse momen- 
tum are mostly produced in hard scattering processes and 
modified by interactions with the medium in high energy 
heavy- ion collisions [72|, [zl, [74| , the change of (j> spectra 
from the Levy function shape in peripheral Au+Au colli- 
sions to an exponential function shape in central Au+Au 
collisions may indicate that different physics dominates 
the particle production in this px region. In the low px 
(pr < 1-0 GeV/c) region, the N part /2 normalized <j> yield 
in d+Au collisions scales with that in p + p collisions, 
whereas it is enhanced significantly in central Au+Au 
collisions. This indicates that the hot environment cre- 
ated by central Au+Au collisions favors the production 
of soft (j> mesons. 

Theoretical calculations have shown that particles with 
different transverse momenta (or in different collision sys- 
tems) are produced by or evolve with different mecha- 
nisms, such as hydrodynamics l35l. [361. [37L HH, coales- 
cence/recombination [33. |40L fall |42|. 43i,E|, fragmenta- 
tion [75|, [zl, [z3] and jet quenching 72|, 173 . [74[ mecha- 
nisms. The observed change of the <fi px spectra shape in 
our measurements is likely due to the change of these pro- 
duction mechanisms in different kinematic regions and 
collision systems. Further evidence of this will be dis- 
cussed later, based on the measurements of different ob- 
servables. 

Figure [T2] presents the 4> meson midrapidity yield per 
participant pair (dN/dy)/(0.5N part ) as a function of 
N P art (approximately proportional to the size of the col- 
lision system). The measured midrapidity yield per par- 
ticipant pair increases nonlinear ly with N partl except for 
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FIG. 10: (Color online) <f> meson transverse mass distributions for different collision systems and different energies. For clarity, distributions 
for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) 
function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the <j> meson spectra for 
Au+Au collisions at 200 GeV in Ref. [50ll to correct for the kaon identification efficiency effect missed previously. 



the largest centrality bins and the Au+Au 130 GeV re- 
sults where there are only three centrality bins with big 
error bars due to the limited statistics. For 200 GeV 
collisions, the yield increases rapidly from p + p and 
d+Au to peripheral Au+Au collisions and then saturates 
for midcentral Au+Au collisions. For the same N part , 
(dN/dy)/(0.5N part ) increases with the collision energy 
of the Au+Au collisions. This is expected because of 
the increase of energy available to produce the <fi mesons. 
The centrality and energy dependences of the enhance- 



ment of (j) meson production can reflect the mechanism 
of strangeness enhancement in a dense medium formed 
in high energy heavy- ion collisions p78j . 

The upper panel of Fig. fTBI shows the N part dependence 
of <j> meson (px) in different collision systems, where (pt) 
is extracted from the best fit to the spectra as de- 
scribed and shown in Table III. The measured (pr) of 
the 4> meson shows no significant centrality dependence 
within systematic errors. At the same N part value, the <fi 
meson (pt) increases slightly with collision energy from 
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62.4 to 200 GeV. 

The mean values of transverse momentum (p?) as a 
function of hadron mass from 62.4 and 200 GeV cen- 
tral Au+Au collisions are presented in the lower panel of 
Fig. [131 These data are taken from Refs. [H EJ. The 
(pr) of ordinary hadrons tt~ , K~ , and p follows a trend 
that is increasing with the mass of the hadron, as ex- 
pected from the dynamics of these particles coming from 
a common radial velocity field shown as the hatched band 
[Hydro, model (1)] in the plot [IE HI]]- However, heavy 
hypcrons such as H and Q, show a deviation from the 
trend. Their values of (pr) are lower than the expected 
ones. The observed (px) values for (j) meson and A are 
similar to those of 3 and £1. Meanwhile, another hydro- 
dynamic model (2) shown by the curve [IE EH , which 
considers possible different chemical freeze-out temper- 
atures for ordinary and strange hadrons, gives a better 
description for strange particle (pt)- This behavior can 
be explained if strange hadrons have a smaller scattering 
cross section than ordinary hadrons in the later hadronic 
stage of the collisions. These strange particles would 
then decouple earlier from the system. The collective 
motion of the <f> meson and multistrange hadrons S and 
f2 should have been developed at the early partonic stage 
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FIG. 12: (Color online) Npart dependence of (dN/dy)/ (0.5N par t) 
in five different collision systems: Au+Au 62.4, f30, 200 GeV; p+p 
200 GeV (Inelastic); and <i+Au 200 GeV. Statistical and systematic 
errors arc included. 



in Au+Au collisions at RHIC. If radial flow is built up 
through the evolution of the system, the particles with a 
smaller hadronic cross section would have smaller radial 
velocity and relatively smaller (px). 
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FIG. 11: (Color online) Comparison of transverse momentum 
spectra shape among different 200 GeV collision systems: Au+Au 
(0-5%), d+Au (0-20%) and p + p (inelastic)). The spectra are nor- 
malized by N bin (top panel) and N par t/2 (bottom panel). 



C. Ratios 

The yield ratio 4>/tt~ as a function of the center-of- 
mass energy per nucleon pair (y/s NN ) is presented in the 
upper panel of Fig. [Til The <t>/ir~ ratio increases with 
energy in both A+ A [ij, HE El HE [H| and p+p col- 
lisions (13, HE HI, Ia3| > indicating that the yield of the cf) 
increases faster than that of the ir~ in A+A and p+p col- 
lisions with increasing ^/s NN . The 4>/tt~ ratio in A+A 
collisions is enhanced compared to that for p + p colli- 
sions, which indicates that A+A collisions may provide 
a more advantageous environment for the production of 
<f> mesons. In fact, an enhanced production of (f> meson 
in heavy-ion environment has been predicted to be a sig- 
nal of QGP formation [f|. However, no clear conclusion 
can be drawn from the experimental measurements since 
the relative enhancement from p + p to A+A collisions 
does not seem to change for y/s^ >10 GeV. The bot- 
tom panel of Fig. [F|] shows the N part dependence of the 
4>/tt~ ratio in different collisions. The 4>/t:~ ratio first 
increases with N part , and then seems to be saturated in 
the high N part region. 

To further study whether (j) meson production, or just 
that of strange particles, is enhanced in high energy A+A 
collisions with respect to elementary collisions, we have 
plotted the yield ratio of 4>/K~ as a function of ^/s^^ in 
A+A, e + e and p+p collisions in the top panel of Fig.[T5l 
For these collisions, at energies above the threshold for 
<f> production, the 4>/K~ ratio is essentially independent 
of collisio n sp ecies and e nergy f rom a few GeV up to 

200 Gev [3 [H EH, H, HI IE H3, HE HE Hi- The 
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TABLE III: Results from fits to the transverse mass distributions of the <f> meson. The fit functions used to extract the results 
are also listed. All values are for midrapidity \y\ <0.5. The first error is statistical; the second is systematic. 





Centrality 


Fit Function 


X 2 /ndf 


T exp/Levy (McV) 


n 


<p T > (MeV/c) 


dN/dy 


Au+Au 


0-20% 


Exp. 


8.4/9 


328± 6 ± 22 


- 


922 ± 13 ± 61 


3.52 ± 0.08 ± 0.45 


(62.4 GoV) 


20-40% 


Exp. 


8.4/9 


324± 6 ± 23 


- 


913 ± 12 ± 65 


1.59 ± 0.03 ± 0.15 




40-60% 


Exp. 


14.5/9 


308± 8 ± 25 


- 


881 ± 16 ± 71 


0.58 ± 0.01 ± 0.07 




60-80% 


Exp. 


13.3/9 


279± 9 ± 28 


- 


822 ± 19 ± 82 


0.15 ± 0.004 ± 0.02 


Au+Au 


0-11% 


Exp. 


5.3/7 


379+ 50 ± 45 


- 


1095 ± 147 ± 131 


5.73 ± 0.37 ± 0.57 


(130 GoV [49]) 


11-26% 


Exp. 


3.2/5 


369+ 73 ± 44 


- 


1001 ± 144 ± 120 


3.33 ± 0.38 ± 0.33 




26-85% 


Exp. 


9.0/6 


417+ 75 ± 50 


- 


1021 ± 99 ± 123 


0.98 ± 0.12 ± 0.10 


Au+Au 


0-5% 


Exp. 


11.0/12 


357 ± 3 ± 23 


- 


977 ± 7 ± 64 


7.95 ± 0.11 ± 0.73 


(200 GcV) 


0-10% 


Exp. 


10.2/12 


359 ± 5 ± 24 


- 


979 ± 20 ± 66 


7.42 ± 0.14 ± 0.68 




10-20% 


Exp. 


9.7/12 


373 ± 4 ± 26 


- 


1010 ± 8 ± 69 


5.37 ± 0.09 ± 0.50 




20-30% 


Exp. 


26.7/12 


387 ± 4 ± 26 


- 


1022 ± 14 ± 68 


3.47 ± 0.06 ± 0.44 




30-40% 


Exp. 


21.1/12 


371 ± 4 ± 24 


- 


1005 ± 8 ± 64 


2.29 ± 0.04 ± 0.23 




40-50% 


Levy 


17.4/11 


315 ±11 + 38 


22.7 ± 4.3 


949 ± 13 ± 67 


1.44 ± 0.03 ± 0.14 




50 60% 




6 9/11 


9qn 4- n -4- 


13 8 + 1 9 




n 89 4- n 09 -1- 09 




60-70% 


Levy 


7.4/11 


291 ±13 + 29 


18.6 ± 3.6 


926 ± 15 ± 75 


0.45 ± 0.01 ± 0.05 




70-80% 


Levy 


5.5/11 


243 ± 15 ± 25 


13.0 ± 2.3 


851 ± 19 ± 85 


0.20 ± 0.01 ± 0.02 


p + p (200 GcV, NSD [14]) 


0-100% 


Levy 


10.1/10 


202 ± 14 ± 11 


8.3 ± 1.2 


812 ± 30 ± 41 


0.018 ± 0.001 ± 0.003 


d+Au 


0-20% 


Levy 


4.7/11 


323 ± 20 ± 32 


15.5 ± 3.9 


1030 ± 57 ± 103 


0.146 ± 0.005 ± 0.014 


(200 GcV) 


20-40% 


Levy 


13.6/11 


316 ± 19 ± 32 


16.9 ± 4.7 


1007 ± 35 ± 101 


0.103 ± 0.003 ± 0.010 




40-100% 


Levy 


12.4/11 


263 ± 15 ± 26 


12.2 ± 2.1 


920 ± 35 ± 92 


0.040 ± 0.001 ± 0.004 




0-100% (MB) 


Levy 


17.5/11 


297 ±11 + 30 


13.9 ± 1.8 


973 ± 26 ± 97 


0.071 ± 0.001 ± 0.005 



lower panel of Fig. [TBI shows that the yield ratio <fi/K~ 
from our analysis is also almost constant as a function 
of centrality. This is remarkable considering that the 
environment created by p + p collisions is so drastically 
different from that of Au+Au collisions that have both 
partonic and hadronic interactions. 

The centrality dependence of the 4>/K~ ratio provides 
another serious test for rescattering models based on the 
assumption that kaon coalescence is the dominant mech- 
anism for 4> production. The KK and if-hyperon modes 
are included in these rescattering models [±2|, [H| ■ They 
predict an increasing ratio vs. centrality (shown 

by the dashed line in the lower panel of Fig. [T5| , which is 
again in contradiction with the approximately flat trend 
of our measurements. Note also that the <f)/K~ ratio 
from the UrQMD model is scaled by a factor of 3.4 to 
match the magnitude of our measurement for peripheral 
Au+Au collisions. The comparisons of the data to pre- 
dictions from these rescattering models including (pt) 
and <fi/K~ effectively rule out kaon coalescence as the 
dominant production mechanism for the <p meson. These 
measurements of the 4>/ 1 K~ ratio may point to a com- 
mon underlying production mechanism for <p and strange 
mesons in all collision systems. 

In addition, statistical models [HI, H H ^ 
based on the assumption that the accessible phase space 
is fully saturated and thermalized can reproduce the 
STAR measurements of integrated hadron yield ratios 
(including <fi/iv~ and 4>/ K~) at midrapidity in p + p and 
Au+Au collisions at 200 GcV [34]. In these statisti- 
cal models the strangeness phase-space occupancy factor 
7 S approaches unity with increasing N part , which indi- 



cates that strangeness approaches/reaches chemical equi- 
libration for midcentral and central Au+Au collisions at 
RHIC energies. 

Because the mechanisms of (multi)strangc particle pro- 
duction are predicted to be very sensitive to the early 
phase of nuclear collisions, the ratio of f2/</> is expected 
to reflect the partonic nature of the thermal source that 
characterizesQGP [96[ and the effects of the strong color 
field (SCF) [98|. In Fig. [lj the ratios of Q/(j> vs. p T are 

presented for different centrality bins. The fl + 
data poi nts arc from Ref. [3J] (for 0-12%) and from 
ref. [lOOj for the other centralities. Also shown in the 
figure are three curves from two recombination model 
expectations for central collisions. A model by Ko et al., 
based on the dynamical recombination of quarks [97j , is 
compared with the data and found to overpredict the 
ratio by a factor of about 3 over the whole pt region 
(Note: the bumpy shape for this model is due to the 
limited statistics of their Monte Carlo simulation of the 
model.). Based on <f> and f2 production from coalescence 
of thermal s quarks in the medium (9|| , Hwa and Yang 
can describe the trend of the data up to pt ~ 4 GeV/c 
but fail at higher px (solid line) . In the alternative SCF 
scenario [H, a large string tension of k = 3 GeV/fm 
can reproduce the data up to pt ~ 4.5 GeV/c in the 
framework of the HIJING/BB v2.0 model, but the ef- 
fect of strong color electric fields remains an open issue 
and needs to be further investigated. With decreasing 
centrality, the observed f2/</> ratios seem to turn over 
at successively lower values of p<r, possibly indicating a 
smaller contribution from thermal quark coalescence in 
more peripheral collisions (9o| . This is also reflected in 
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FIG. 13: (Color online) Top panel: N par t dependence of (pt)<p in 
different collision systems; Bottom panel: Hadron mass dependence 
of (pt) in central Au+Au collisions at 62.4 and 200 GeV. The band 
and curve show two hydrodynamic model calculations for central 
Au+Au collisions at 200 GeV. Note: Hadron masses for the Au+Au 
62.4 GeV data are shifted slightly in the x-axis direction for clarity, 
and systematic errors are included for the <j>. 



the smooth evolution of the pr spectra shapes from the 
thermal-like exponentials to Levy-like curves. 

D. Nuclear modification factor 

The measurement of the nuclear modification factors 
R cp and Rab provides a sensitive tool to probe the pro- 
duction dynamics and hadronization process in relativis- 
ts heavy-ion collisions [ii mm mm si- Rcp ' which 

is the ratio of yields in central to peripheral heavy-ion 
collisions normalized by A^„, is defined as 



RcpiPT) 



[dN/(N bin dp T )} 



central 



[dN / (N bin dp T )} peripheral ' 



(14) 



and Rab j which is the yield ratio of nucleus (A) + nu- 
cleus (B) collisions to inelastic p+p collisions normalized 
by Nbin, is defined as 



Rab(pi 



[dN/(N bin dp T )} 



A+B 



[dN/dp 7 



}p+p 



(15) 



FIG. 14: (Color online) Top panel: Energy dependence of the ratio 
(f>/n~ in A+A (full points) and p + p (open points) collisions. Stars 
are data from the STAR experiment at RHIC. Bottom panel: N par t 
dependence of ratio 4>/tt~ in different collision systems. Systematic 
errors are included for the STAR data points. 



where Ntn„ is the number of binary inelastic nucleon- 
nucleon col lision s determined from Glauber model calcu- 
lations [xJ, llOlj . It is obvious that these two ratios will 
be unity if nucleus-nucleus collisions are just simple su- 
perpositions of nucleon-nucleon collisions. Deviation of 
these ratios from unity would imply contributions from 
nuclear or QGP effects. It should be mentioned that 
when using p + p collisions in the Rab ratio, inelastic 
collisions should be used instead of the measured NSD 
collisions, therefore a correction factor 30/42, which was 
discussed in the above subsection B, was applied to cor- 
rect the NSD yield to the inelastic yield in p+p collisions. 
There is a py-dependent correction to the NSD distribu- 
tion from singly diffractive (SD) events, which is small in 
our measured pr range, being 1.05 at pr = 0.4 GeV/c 
and unity above 1.2 GeV/c as determined from PYTHIA 
simulations (7l| . 

Figure [T7| presents the px dependence of R cp for the 
4> with respect to midperipheral (top panel) and most- 
peripheral (bottom panel) bins in Au+Au collisions at 
200 GeV. Results for the A + A, the K° s (H, the tt+ +ir- 
and the p + p [l02j particles are also shown in the figure 
for comparison. Both statistical and systematic errors 
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FIG. 15: (Color online) Top panel: Energy dependence of ratio 
4>/K~ in A+A (full symbols) and elementary (e + e: open triangles 
and p + p: open circles and star) collisions. Stars are data from 
STAR experiments at RHIC. Bottom panel: N par t dependence of 
ratio 4>/K~ in different collision systems. The dashed line shows 
results from UrQMD model calculations. Systematic errors are 
included for the STAR data points. 



are included. Most of the systematic errors cancel in the 
ratios; however, the uncertainty due to particle identifi- 
cation from dE/dx remains as the dominant source and 
varies from point to point over the range ~ 7% — 12%. 
In the measured px region, the R cp of <f> meson is con- 
sistent with Npart scaling at lowest px (dot-dashed line) , 
and is significantly suppressed relative to the binary col- 
lision scale (dashed horizontal line at unity) at all pt in 
Au+Au collisions at 200 GeV. When compared to the 



STAR measured A + A and K s data [4f 



the R cp for 



follows that of Kg rather than that of the similarly mas- 
sive A (A), especially for the case of 0-5%/40-60%. For 
0-5%/60-80%, the R cp of sits between that for the K s 
and the A. This may be attributed to the shape change 
of the 4> spectra from exponential at 40-60% centrality to 
Levy at 60-80% centrality as discussed above, which may 
be due to the change of the </> production mechanism at 
intermediate pt in different environments with different 
degrees of strangeness equilibration. 

Figure [TS] presents the pt dependence of R cp for the 
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FIG. 16: (Color online) The Q/<f> ratio vs. px for three centrality 
bins in ^/sjvjv = 200 GeV Au+Au collisions, where the data points 
for 40-60% are shifted slightly for clarity. As shown in the legend, 
the lines represent the results from Hwa and Yang [96| , Ko et al. [97j 
and for SCF, Refs. iHHl. 



4> meson for Au+Au 62.4, 200 GeV and d+Au 200 GeV 
collisions. For the three collision systems, the R cp factor 
follows N par t scaling at low px- However, the R cp at in- 
termediate pt is strongly suppressed in Au+Au collisions 
at 200 GeV, is weakly suppressed in Au+Au collisions at 
62.4 GeV, and shows no suppression in <i+Au collisions 
at 200 GeV. 

Our measured R cp for the (j> meson further supports the 
proposed partonic coalescence/ recombination scenario at 
intermediate px 0, IH, S3, |48[ , where the centrality de- 
pendence of particle yield depends on the number of con- 
stituent quarks (NCQ) rather than on the mass of the 
particle. Das and Hwa [47j proposed a recombination 
scenario to explain the leading particle effect in p + p col- 
lisions, where the fragmenting partons would recombinc 
with sea quarks to form hadrons. The resulting px dis- 
tribution for the leading particle is therefore determined 
by the fragmenting quarks. In heavy-ion collisions, it is 
possible that qq (qqq) come together and form a meson 
(baryon) due to the high density of quarks (antiquarks) in 
the collision system. The abundant nearby partons give 
the recombination/coalescence mechanism a comparative 
advantage over the fragmentation mechanism for particle 
production in the intermediate px region, and this suc- 
cessfully explains the relative enhancement of baryons in 
the intermediate Px r egion, such as the ratios p/w and 



A/K s [lOl [lol Il04|. Recently, Hwa and Yang used 
their recombination model to successfully describe the 4> 
meson px spectra, but it fails to reproduce the px de- 
pendence of fl/4> at higher px as shown in the previ- 
ous subsection. They assumed that for the production 
of <p this thermal component dominates over others in- 
volving jet shower contributions because of the suppres- 
sion of shower s quarks in central Au+Au collisions [96| . 
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(Color online) py dependence of the nuclear modification 
i cp in Au+Au 200 GeV collisions. The top and bottom pan- 



FIG. 17: 

factor R c 

els present R cp from midperipheral and most-peripheral collisions, 
respectively. See legend for symbol and line designations. The 
rectangular bands show the uncertainties of binary and participant 
scalings. Statistical and systematic errors are included. 



However, the partonic medium, with an intrinsic tem- 
perature, may modify the recombination probability and 
momentum distribution for the <fi }42l . l48j . 

Figure Q33 presents the px dependence of the nuclear 
modification factor Rab in Au+Au and d+Au collisions 
at 200 GeV. For comparison, data points for RdAu of 
7r + + 7r~ and p + p are also shown in the figure. The 
RdAu of 4> mesons reveals a similar enhancement trend as 
those of 7r + +7r~ and p + p at the interme diat e pt , which 
was attributed to be the Cronin effect KM, [iM H07l. 



The Cronin enhancement may result ei ther from mo- 
mentum broade ning due to multiple soft [l08| (or semi- 
hard pH, [Tl3 HH [113) scattering in the initial state 
or from final st ate in teractions suggested in the recom- 
bination model [l!3j |. These mechanisms lead to differ- 
ent particle type and/or mass dependence in the nuclear 
modification factors as a function of pt- Our measure- 
ment of RdAu of 4> mesons does not have th e precision 
to differentiate particle dependence scenarios |ll4 , Ill5l | . 
On the other hand, the Raa in Au+Au 200 GeV is lower 
than that in d+Au 200 GeV and consistent with the bi- 
nary collision scaling at intermediate px- These features 
are consistent with the scenario that features the onset 



CC 



10" 




Scaling — binary 



participant 



Au+Au 200 GeV Au+Au 62.4 GeV d+Au 200 GeV 



0-5% 
60-80% ' 



0-20% 
60-80% 



0-20% 
40-100% 



□ 



1 i 1 1 1 1 i 
4 5 

p T (GeV/c) 



FIG. 18: (Color online) px dependence of the nuclear modifica- 
tion factor R cp in Au+Au 62.4 and 200 GeV and d+Au 200 GeV 
collisions, where rectangular bands represent the uncertainties of 
binary and participant scalings (see legend). Statistical and sys- 
tematic errors are included. 
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FIG. 19: (Color online) pt dependence of the nuclear modification 
factor R AB for <j> in Au+Au 200 GeV and d+Au 200 GeV collisions. 
For comparison, data points for 7r++7r _ in d+Au 200 GeV and p+p 
in d+Au 200 GeV are also shown (see legend). Rectangular bands 
shows the uncertainties of binary (solid line) and participant (dot- 
dash line) scalings. Systematic errors are included for if>, 7r + + tt~ 
and p + p. 



of parton-medium final state interactions in Au+Au col- 
lisions. The two Rab observations in central d+Au and 
Au+Au collisions are consistent with the shape compar- 
isons in Fig. [TT] (top panel). 

E. Elliptic flow 

Figure [50] shows the elliptic flow of the cf> meson 
as a function of pt in MB (0-80%) Au+Au collisions at 
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FIG. 20: (Color online) px dependence of the elliptic flow V2 of <j>, 
A, and in Au+Au collisions (0-80%) at 200 GeV. Data points 
for <j> are from the reaction plane method (full up-triangles) and 
invariant mass method (full circles), where data points from the 
reaction plane method are shifted slightly along the x axis for clar- 
ity. Vertical error bars represent statistical errors, while the square 
bands represent systematic uncertainties. The magenta curved 
band represents the V2 of the 4> m eson from the AMPT model with 
a string melting mechanism Ill6ll . The dash and dot curves rep- 
resent par amet rizations inspired by number-of-quark scaling ideas 
from Ref. [TPl for NQ=2 and NQ=3 respectively. 



^s NN =200 GeV. The w 2 of the K s and A measured by 
STAR in Au+Au 200 GeV collisions are also shown 
for comparison. Measurements of V2 for the (ft meson 
from both reaction plane and invariant mass methods 
are presented, and they are consistent with each other. 

The first interesting observation is that the (ft meson 
has significantly nonzero 1)2 in the measured px region. 
If the (ft meson has a small interaction cross section with 
the evolving hot-dense matter in A+A collisions, it would 
not participate in the late-stage hadronic interactions in 
contrast to hadrons such as tt, K, and p(p) which freeze- 
out later. This indicates that the nonzero V2 of the (ft 
meson must have been developed in the earlier partonic 
stage. In the low px region (<2 GeV/c), the v 2 value 
of (ft is between that for the K s and the A in Au+Au 
200 GeV collisions, consistent with the expectation of a 
mass ordering for vi in hydrodynamic models. These 
observations support the hypothesis of the development 
of partonic collectivity and possible thermalization i n th e 
early stages of heavy- ion collisions at RHIC [H, 1 1 OOj ] . 
although the underlying mechanism for the equilibration 
process remains an open issue. 

In the intermediate pr region (~2-5 GeV/c), the V2 of 
the (ft meson is consistent with that for the K s rather 
than for the A. When we fit the t^ far) of <ft mesons with 
the quark number scaling function [117| , the resulting fit 
parameter NCQ (number of constitute quarks) = 2.3 ± 
0.4. The fact that the vi(j>t) of (ft is the same as that 
of other mesons indicates that the heavier s quarks flow 



FIG. 21: (Color online) Elliptic flow 1)2 as a function of px, ^2(pt), 
for the tj> meson from different centralities. The vertical error bars 
represent the statistical errors while the square bands represent 
the systematic uncertainties. For clarity, data points of 10-40% are 
shifted in the px direction slightly. 

as strongly as the lighter u and d quarks. The AMPT 
model with string melting and parton coalescence mech- 
anisms can reproduce the experimental results well up 
to 3 GeV/c, which favors the hadro nizat i on s cenario of 
coalescence/recombination of quarks |ll6tlll8j . 

The V2 of the (ft meson from other centralities arc shown 
in Fig. 1211 The data are analyzed from the invariant mass 
method only. As expected, v 2 {pt) increases with increas- 
ing eccentricity (decreasing centrality) of the initial over- 
lap region. This trend is also illustrated in Table IIVI 
which presents the pT-integrated values of </>-meson ellip- 
tic flow, (1)2) , calculated by convoluting the V2{pt) with 
the respective px spectrum for four centrality bins. It 
should be noted that the centrality dependence of the 
(^2) of the (ft meson is consistent with that of the charged 
hadrons [671 ] . 



TABLE IV: Integrated elliptic flow (1J2) for the (ft meson for 
four centrality bins in Au+Au collisions at 200 GeV. 



Centrality (%) 


{v 2 } (%) 


40 - 80 


8.5 ± 1.1 (stat) ± 0.2 (sys) 


10 - 40 


6.6 ± 0.8 (stat) ± 0.2 (sys) 


0-5 


2.1 ± 1.2 (stat) ± 0.5 (sys) 


0-80 


5.8 ± 0.6 (stat) ± 0.2 (sys) 



IV. CONCLUSION 

In conclusion, STAR has measured (ft meson produc- 
tion for 62.4, 130, 200 GeV Au+Au, 200 GeV d+Au, 
and NSD p + p collisions at RHIC. Details of the analysis 
method for (ft meson are presented. The respective en- 
ergy and N part dependence of the (ft meson production, 
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as well as the px spectra for five collision systems, are 
reported. 

The 4> spectra in central Au+Au 200 GcV collisions 
are described well by an exponential function. The spec- 
tra for p + p, d+Au and the most peripheral Au+Au 
200 GcV collisions are better described by a Levy func- 
tion due to the high-p^ power-law tails. This change of 
spectra shape from p + p, d+Au and peripheral Au+Au 
to central Au+Au collisions is most likely due to a change 
of the dominant <p production mechanism in the different 
collision environments. The yield of <j> mesons per par- 
ticipant pair increases and saturates with the increase of 
N P art- It is found that the (px) of <fi, A, 5, and f2 does 
not follow the (pt) vs. hadron mass trend determined 
by the 7r~, K~ , and p. This may be due to their small 
hadronic cross sections, which indicates that the </> and 
strange hadrons can retain more information about the 
early state of the collision system. 

The 4>/K~ yield ratios from p + p and A+B colli- 
sions over a broad range of collision energies above the 
<f> threshold are remarkably close to each other, indicat- 
ing similar underlying hadronization processes for soft 
strange quark pairs in these collisions. The lack of a 
significant centrality dependence of the <£>/ 1 K~ yield ra- 
tio and (pr)<t> effectively rules out kaon coalescence as a 
dominant production channel at RHIC. The trend of the 
£l/<j> ratio is consistent with recombination models and a 
strong color field scenario up to px ~ 4 GeV/c in central 
Au+Au collisions. 

The measurement of the <p meson nuclear modification 
factor Rab is consistent with the Cronin effect in d+Au 
200 GcV collisions, and with the energy loss mechanism 



in Au+Au 200 GeV collisions. The <f> meson was found 
to have nonzero vi in the measured px range. When 
comparing the <f> meson nuclear modification factor R cp 
(0-5%/40-60%) and elliptic flow to those of the simi- 
lar mass A baryon, and to the lighter Kg meson, we see 
that the <j> meson clearly behaves more like the K s me- 
son than the A baryon. Therefore, the processes relevant 
to R cp and vi at intermediate px are driven not by the 
mass of the particle, but rather by the type of the par- 
ticle, i.e. number of constituent quarks (NCQ) scaling. 
The coalescence/recombination model provides a fairly 
consistent picture to describe particle production in the 
intermediate px region over a broad range of collision 
energies and system sizes at RHIC. 
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